CURR£NT R ~Jrpfn 4> CRGSS ~TAIL 
Sfy> k FIGL °- A1 -ISNeO CURREM 
(Nev« Hampshire Univ*) < 


AND 

P 


N95-19032 
Unci as 


£3/46 0038395 





Abstract. Orbits of individual charged particles were traced in a one-dimensional magnetic 
field model that included a uniform cross-tail component B yo . The effects of B yo on the cross-tail 
current distribution j y (z ) , the average cross-tail drift velocity (v (z ) ), and the average pitch 
angle change (Aa) experienced during current sheet encounters were calculated. The addition of 
a B yo that exceeded several tenths of one nanotesla completely eliminated all resonance effects 
for odd-A orbits. An od d-N resonance involves ions that enter and exit the current sheet on the 
same side. Pitch angles of nearly all such ions changed substantially during a typical current sheet 
interaction, and there was no region of large cross-tail drift velocity in the presence of a modest 
B yo . The addition of a very large B yo guide field in the direction that enhances the natural drift 
produces a large j y and small (Aoc) for ions with all energies. The addition of a modest B yo had 
less effect near even-A resonances. In this case, ions in a small energy range were found to 
undergo so little change in pitch angle that particles which originated in the ionosphere would 
pass through the current sheet and return to the conjugate ionosphere. Finally, the cross-tail drift 
of ions from regions dominated by stochastic orbits to regions dominated by either resonant or 
guiding center orbits was considered. The ion drift speed changed substantially during such tran- 
sitions. The accompanying electrons obey the guiding center equations, so electron drift is more 
uniform. Any difference between gradients in the fluxes associated with electron and ion drifts 
requires the presence of a Birkeland current in order to maintain charge neutrality. This plasma 
sheet region therefore serves as a current generator. The analysis predicts that the resulting Birke- 
land current connects to the lowest altitude equatorial regions in which ions drift to or from a 
point at which stochastic orbits predominate. The proposed mechanism appears only in analyses 
that include non-guiding-center effects. 
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1. Orbit Types and Magnetotail Regions 

A number of studies have traced charged particle orbits in taillike magnetic field models. 
Chen [1992] reviewed this topic with particular emphasis on the dynamics of particle orbits. The 
one-dimensional modified Harris [1962] model 

B x (z) = B xo tmh (z/L) 

B y (z)=0 (J) 

B rW = B zo 

has been employed in many of these analyses. The constant parameters B xo , B zo , and L can be 
selected so that, within the current sheet, the shapes of Harris-type field lines are similar to the 
shapes of magnetotail field lines in the more comprehensive Tsyganenko [1989] or T89 model 
[Kaufmann et al., 1993 a, hereafter referred to as paper 1]. Since one-dimensional models cannot 
have a point at which B z = 0, they are not able to represent the midtail plasma sheet boundary 
layer (PSBL) or lobes. In this work, we investigate particle motion only within the inner portion 
of the central plasma sheet (CPS). Equation (1) was used at Izl < L and 


B x( z ) 


B, 




tanh [ (z/L) ^] + 5 


1 -- 


tanh ( 1 ) 


( 2 ) 


with rj = 5 was used when Izl > L. Equation (2) was selected because it provides a smooth cutoff in 
the cross-tail electric current density, j y (z) = dB x Idz , between Izl = L and Izl = 4 L/ 3 . Beyond 

this point, which generally is well within the plasma sheet, particles carry little cross-tail current 
and orbit tracing can be stopped. 

Invariant Tori 


A charged particle moves along the surface of a torus in phase space, or of a torus in ordinary 
x, y, z coordinates for the one-dimensional magnetic field considered here [Buchner and Zelenyi , 
1989]. The entire torus is ring-like for a class of ions that remain trapped within the principal cur- 
rent sheet (Figure la). These ions continually bounce or meander back and forth across z = 0. 
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When projected onto the y-z plane, the ring-like or crossing trajectories resemble a figure 8. This 
orbital shape is especially important in studies of j y ( z ) . Ring-like or figure 8 ions carry current in 
the positive y direction at the largest Izl they reach, and in the negative y direction near z = 0. 

For some orbits the earthward part of a torus (right side of Figures lb and lc) splits into two 
segments, one in the northern and one in the southern hemisphere. The curve on which the torus 
splits is called an orbital separatrix. An ion bounces back and forth across z — 0 while on the ring- 
like segment (left sides of Figures lb and lc) and spirals around a magnetic field line while on the 
split segments of a torus. This ion crosses one separatrix when it enters the ring-like segment and 
one when it leaves. Depending on its pitch angle, the ion either may mirror during the spiraling 
portion of the orbit, or it may escape in a one-dimensional model. Mirroring trajectories also are 
referred to as cucumber orbits. Such ions are trapped in the sense that they remain within the CPS, 
but differ from the figure 8 or ring-type trapped orbits because mirroring ions do not remain 
within the principal current sheet. 

Orbital Parameters 

The magnetic moment , | i = mv 2 J (2 B) is commonly used as an approximate invariant of the 
motion when nonrelativistic ions spiral around magnetic field lines. Such particles are character- 
ized by Q » co^, where Q, is the local gyrofrequency and co^ is the bounce frequency between mir- 
ror points. The guiding center approximations are nearly obeyed at large Izl on the split or spiral 
segments of a torus, so that p is nearly conserved on this non-crossing portion of an orbit. 

The action integral I z = (2k) P z dz [Speiser, 1970; Sonnerup, 1971] is approximately 
conserved while particles bounce back and forth across z = 0 on the ring-like or crossing portion 
of an orbit. This part of an orbit is characterized by cd z < Q zo , where co z is the north-south bounce 
frequency across the z = 0 plane and Q zo = qB z Jm is the gyrofrequency based on B zo . In gen- 
eral, neither / z nor p is conserved during separatrix crossings. 

Buchner and Zelenyi [1986; 1989] introduced the K parameter 
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which, together with l z or j. i , helps to organize the dominant trajectories found in a given region of 
the tail. In the model field used here, k 2 = RJ where R c is the magnetic field line radius of cur- 


vature at z - 0 and p z0 - v/Q zq . Thermal ions and electrons have K > 2 in the radiation belts and 


near-Earth tail. These particles spiral around magnetic field lines and obey the guiding center 
equations during all or almost all of an orbit. There is at most a very short ring-like region near z = 
0 where the northern and southern segments join. Particles with K > 2 only cross the z = 0 plane 

once during each current sheet interaction, and they drift from dawn to dusk during this brief 
period. 


Guiding center motion breaks down for most ions beyond a point that usually is located some- 
where between 5 R E and 15 R E in the magnetotail [Pulkkinen et al., 1992; paper 1]. We refer to 
that portion of the non-guiding-center region in which 2 > k > K r = B z JB xo as the midtail. It is 
convenient to divide the associated ions into categories according to k. Groups of ions with differ- 
ent k carry substantially different cross-tail current distributions j (z) , and suffer different aver- 
age pitch angle changes (Aa) during a current sheet interaction. The primary goal of this work is 
to study j y (k, z) and (Aa (k) ). These two parameters were selected because they can be evalu- 
ated using satellite data, and because of their importance in determining the structure of the 
plasma sheet. 

Stochastic and Resonant Orbits 


Both p and l z of nearly all ions with certain values of k (e.g., k = 0.82, 0.40, and 0.27) are sub- 
stantially different before and after almost every current sheet interaction. These orbits are called 
stochastic or chaotic. Most stochastic ions eventually leave the current sheet with such a small 
pitch angle that they escape from a one-dimensional model field. The adiabatic invariants p and I z 
change primarily during the pair of separatrix crossings that occur as particles move from a split 
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torus to a ring-like orbit segment and then back. 

A second class of ions which follow split tori in the midtail are called resonant. The magnetic 
moment and I z of these ions also change substantially during most orbital separatrix crossings. 
However, at certain resonant values of K (e.g., k = 0.53, 0.32, and 0.23) the changes in jll and I z at 
a pair of crossings nearly cancel for a wide range of incident pitch and phase angles [Buchner, 
1991]. Resonant ions tend to follow very simple, predictable orbits. For example, nearly all k = 
0.53 ions cross the z - 0 plane just twice during each interaction with the current sheet. These ions 
appear to simply bounce off the current sheet. Figure lb was generated by following a single K = 
0.53 ion while it mirrored and returned several times to the outer CPS. It was necessary to use a 
resonant ion to illustrate this structure so that nearly the same torus was followed after each 
encounter with the current sheet. The K = 0.53 ion in Figure lb was injected at z > 0, so it always 
entered and left the current sheet on the northern side. Figure lb therefore illustrates only half of a 
full split torus. An ion injected at z < 0 would follow a similar trajectory that entered and left the 
southern side of the current sheet. Figure 1 c was prepared by simply adding a reflected duplicate 
of the trajectory in Figure lb. Figure lc therefore illustrates a full split torus. Figures lb and 1 c 
are shown from different viewpoints to more clearly illustrate all aspects of the torus. 

Ions injected at the northern edge of the CPS at the next resonance (k = 0.32) pass through z = 

0 three times and then move into the southern outer CPS. Such ions would alternately follow the 
northern and southern segments of a torus similar to that in Figure lc. However, whenever the 
changes in [i and I z at a pair of separatrix crossings do not exactly cancel, the ion will mirror at 
different locations and therefore be on a different torus during successive orbits. We found it is 
difficult to identify k = 0.32 ions for which cancellation was exact enough over a wide enough 
range of pitch and phase angles so that a single torus would be followed throughout several orbits. 

At the Mh resonance the angle variable Q z associated with the action I z changes by A0 z = Nn 
[Biichnei and Zelenyi, 1989; 1990]. With this definition of A, the K > 2 guiding center motion cor- 
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responds to N ~ 0, and k - 0.53 resonant ions have N = 1 . The relationship between N and k is 
approximately 


[Chen, 1992; Ashour-Abdalla et al, 1993; Kaufmann and Lu, 1993c, hereafter referred to as paper 

2 ]. 

Nonresonant Orbits 

The k of a typical thermal ion decreases as one moves farther out the magnetotail. The aver- 
age ion velocity decreases slowly with increasing distance [Baumjohann et al., 1989]. In the T89 
Kp - 4 model, L increases slowly and B xo decreases rapidly as one moves out the tail (paper 1). 
For example, L changes approximately 30% while B xo changes by a factor of 2 between x = -10 
and -20 R e . The point k = K r is reached when p xo = L or LqB X0 / (mv) = 1, where p xo is the ion 
gyroradius in a magnetic field B xo . An ion with k > K r does not have enough energy to go directly 
from z = 0 to Izl = L without being either reflected back toward z = 0 or being deflected onto a spi- 
ral orbit. Note that the condition K = K r depends upon L, v, and B rn but not on B since it is B r 
that causes the reflection of an ion. The ring-like segment of an invariant torus, and therefore the 
separatrix, are confined to Izl < L for ions with K > K r . 

When k < K r , orbits following split tori become insensitive to K so that stochastic and reso- 
nant effects are not obvious. We refer to these as nonresonant orbits. Burkhart and Chen [1991] 
showed that resonant peaks in the average stochastic trapping time disappear abruptly at K = K r . 
Nonresonant ions experience relatively little change in a during a current sheet encounter (paper 
1). On the ring-like segment of a split torus, a nonresonant ion with appropriate pitch and phase 
angles can bounce back and forth between north-south reflection points that are well beyond Izl = 
L. If many such ions were present, they would carry substantial j y well beyond the edge of the cur- 
rent sheet. Therefore, only a small number of ions can exist with these pitch and phase angles in 
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any self-consistent current sheet, in which the particles must carry the correct j y ( z ) distribution 
to generate the magnetic field in which their orbits were traced. 

In papers 1 and 2 we studied the pitch angle changes and current carrying properties of vari- 
ous groups of resonant ions in the B yo = 0 modified Harris field. The present paper is primarily a 
study of the effect of nonzero B yo on j y (z) and <Aa>. It is found that changes in B yo can pro- 
duce some of the same macroscopic effects that previously were associated with changes in k 
when B yo - 0. However, we find that the odd- and even-A orbits behave quite differently. The 
most striking effect is that odd-A resonances vanish when B yo exceeds several tenths of one nan- 
otesla. 

Finally, the field-aligned currents which are expected if k or B yo varies along a drift path are 
discussed. This rather speculative field-aligned-current generation mechanism is included because 
of its potential importance to magnetosphere-ionosphere coupling. 

2. Cross-Tail Magnetic Fields 

Oveiwiew with B yo = 0 

Figure 2 shows how j y (z) , the distribution of cross-tail current carried by a group of monoen- 
ergetic ions, depends on K when B yo = 0. Groups of ions with randomly selected pitch and phase 
angles were injected at z t = L into a model magnetic field (B = 32 nT, B ~0 B = 3 3 nT L 
= 1 R e ). As is common in one-dimensional studies, this analysis was carried out in a reference 
frame with E y = 0. Pitch angles were selected to have the same distribution as would be found in a 
bi-Max wellian plasma with T^/T ± - 10, while phase angles were selected from an isotropic dis- 
tribution. The structure in Figure 2 also can be seen when T^/T^_ = 1, but is more clearly defined 
with the larger ratio. We are interested in distributions with large T u /T ± at Izl = L partly because 
they are required to produce force balance in a self-consistent one-dimensional model (paper 2). 
Another reason for injecting groups with large 7^/7^ is to study particles that reach low- and 
middle-altitude satellites. All ions in any one group were injected with the same energy rather 
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than with energies selected from a Maxwellian distribution so they would all have the same k. 

The region Izl < 4L/3 was divided into 40 boxes, each having a width Az = L/ 15. Since our 
model is symmetric, results at positive and negative z were added to yield plasma parameters in 
each of 20 boxes between Izl = 0 and 4L/3 . The density n(z ) and electric current distribution 
j y (z) were evaluated for each group by keeping track of the time spent and the distance each ion 
moved across the tail while it was in a given box (paper 2). Since the absolute value of j (z) 
depends on particle energy, B xo , and B zo , we plot a normalized average cross-tail ion velocity 
instead of j y (z) . The average velocity (v y (z) > is given by j y (z) / [ qn (z) ] , where q is the ion 
charge. The vertical axis in Figure 2 is {v y )/v Q , where v 0 is the total ion velocity. This normalized 
cross-tail velocity peaks at more than 0.5 near z = 0 at each of the first three resonances. The peaks 
extend to larger Izl as k gets smaller. Figure 2 also shows that ( v y )/v 0 is small at all z for stochas- 
tic ions, such as those with K = 0.40 or 0.82. 

Simple models ofB y 

Satellites frequently observe magnetotail fields with large components in the geocentric solar 
magnetospheric y direction. However, not all such observations can be modeled by a modified 
Harris field (1) that includes a constant nonzero B yo . The sketches in Figure 3 illustrate several 
reasons why satellites observe a substantial B y . The remaining sections of this paper describe 
effects of only one of these sources. 

A dipolar field has a large B y everywhere except near noon and midnight. The apparent drag- 
ging of magnetospheric field lines into the tail by a force in the equatorial plane (Figure 3a) 
reduces B y in the midtail. However the distorted field lines in a tail of varying width retain a posi- 
tive B y at y < 0 and a negative B y at y > 0 in the northern hemisphere, as shown. The direction of 
this B y is reversed in the southern hemisphere and is zero at the neutral sheet. The magnetic field 
distortion shown in Figure 3 a is primarily produced by currents inside the Earth and by the mag- 
netopause and cross-tail current systems. Field-aligned sheet currents also contribute to the distor- 
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tions sketched in Figure 3 a. If the region 1 current system flows in the outer CPS or PSBL in the 
midtail, then field lines in the inner CPS will have a larger \B y \ than do field lines in the lobes. 

Figure 3b shows another reason why satellites observe magnetic fields in the solar magneto- 
spheric y direction. The neutral sheet is shown tipped, with B everywhere perpendicular to the 
sheet. Dipole tilt is known to produce substantial warping of the neutral sheet [Fairfield, 1979; 
Tsyganenko, 1989]. Since B x - 0 at z = 0 in (1), this equation can describe the fields in Figure 3b in 
a tipped neutral sheet coordinate system, with z measured normal to the neutral sheet. The particle 
orbits in this tipped coordinate system are correctly described using (1) with B yo - 0. Warping is 
generated by the same current systems that were discussed above. 

Figure 3c shows the only type of B y that will be investigated here. In this case, B is no longer 
normal to the neutral sheet. A uniform B yo is added throughout the entire current sheet to describe 
this field. The two ends of a field line in Figure 3c reach Earth at different longitudes. When B yo > 
0, the northern ionospheric footpoint of a magnetotail field line is duskward of the southern foot- 
point. A B y that is uniform throughout the tail is generated by currents outside the magnetosphere. 
We place our model current sheets at infinity. This type of magnetotail B y usually is described in 
terms of penetration or leakage of the y component of the interplanetary magnetic field (IMF) into 
the magnetosphere. As noted above, region 1 Birkeland sheet currents can produce a relatively 
uniform B y in the inner CPS, and a different relatively uniform B y in the lobes. 

Observed B y 

There have been a number of studies of B y in the midmagnetotail. Measurements also have 
been made at synchronous orbit [Cowley and Hughes, 1983] and in the distant tail [Sibeck et al, 
1985]. We are most interested in the midtail within 60 R E of Earth, where the T89 magnetosphere 
model can be used to estimate the parameters in (1). Fairfield [1979] removed the average effects 
of tail flaring, aberration due to motion of the Earth around the sun, and neutral sheet warping to 
determine the residual B y . A correlation with the measured IMF B y yielded an average 13% pene- 
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tration for the IMP 6 data set, taken throughout the plasma sheet and lobes between x = -20 and 
-33 R e . Lui [1984] used an IMP 6 data set that was restricted to neutral sheet crossings. He 
obtained an average 50% penetration of the IMF B y . Sergeev [1987] used ISEE 1 data between x = 
-15 and -23 R E . He separated plasma sheet and lobe data, finding an average 60% penetration in 
the plasma sheet and 15% in the lobes. This difference in B y between the two regions was attrib- 
uted to the effect of region 1 sheet currents in the PSBL. 

K °y maz et al t 1994 l carri ed out an extensive study of B y using IMP 8 magnetotail data taken 
between x = -25 and -40 R E . They found an average 26% penetration when I Izl <8 R E and 9% 
when izl >8 R E . These z values were corrected for average warping effects. However, this study 
showed that the penetration of B y is more complex than can be modeled by a constant B yo . The 
measured B y was much stronger in the flanks than near midnight. This observation suggests that 

IMF B y penetration may be closely associated with processes that take place within or near the 
LLBL. 

A uniform positive B yo could be modeled by a uniform current sheet flowing earthward at 
large positive z, either outside the northern magnetopause or in the northern mantle, and a similar 
tailward flowing current sheet in the southern hemisphere. One way to model the observed 
enhanced B y in the flanks would be to have the current sheets at large Izl become stronger as lyl 
increases. Another possible model would have a uniform current sheet in the mantle and also 
reversed currents near midnight in the PSBL (tailward in the north and earthward in the south 
when B y is positive). We previously evaluated j n for the T89 Kp = 5 model [Kaufmann et al . , 
I" 3 b]. This model has no average uniform B y because it was based on a large data set taken both 
when the IMF B y was positive and negative. However, it may be significant that the T89 model 
showed the largest magnitudes of at large Izl in the flanks. 

Kaymaz et al. [1994] also found an average neutral sheet tilt (Figure 3 b) that correlated with 
IMF B y . The average tilt angle was 7°. The neutral sheet was tilted in a direction which increases 
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the need for a nonzero model B yo in ( 1 ). For example, when B y was negative, as shown in Figure 
3b, the neutral sheet was found to tip in the direction opposite from that shown in our sketch. This 
same distortion had been found in the distant tail by Sibeck et al. [1985], However, once again the 
Kaymaz et al. [1994] study shows the magnetotail is substantially more complex than our model. 
They found the tilt angle much larger than the 7° average in the flanks and much smaller than the 
average, or even reversed, near midnight. 

The above observations suggest that some B y usually should be added to (1) in the midtail cur- 
rent sheet. The uniform B yo we use may be appropriate when following a single orbit, in which an 
ion moves only 1 R £ in the y direction. However, a substantially more complex B y field will be 
needed in a three-dimensional study. The average IMF \B y \ is approximately 2 nT, so it appears 
that the average uniform part of | B y \ is on the order of 1 nT in the neutral sheet, and several tenths 
of one nanotesla at larger tzl. Larger fields are required in the flanks, and the observed fluctuations 
suggest that several nanotesla B y fields can be found almost anywhere at times. The following sec- 
tions show how a uniform B yo with such magnitudes affects particle dynamics and the structure 
of the current sheet. 


3. Energization 

Constants of the Motion 

In the original Harris [1962] magnetic field model (B yo = B z0 = 0), the Hamiltonian 
^(P> r ) - c { [P - <7 A (r )] 2 + m 2 c 2 } and both the x and y components of the canonical 
momentum [P x = ymv x , P y = ymv y + qA y , where y = (1 -v 2 /c 2 )~' /2 ] are constants of the 
motion. The Harris vector potential usually is written with only a y component 

= ~ B xo L ln t cosh ( z/L ) ] (5) 

The Poisson brackets [H, PJ, [H,P y ], and [P x ,P y ] all are zero, so that the three constants of 
the motion are independent, or are in involution. 

The addition of a nonzero B zo to an infinite one-dimensional current sheet, as in (1) and (2), 
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can be described by modifying A xt A y , or both. The dynamics of orbits in this modified Harris 
model were studied by Chen and Palmadesso [1986] and by Buchner andZelenyi [1989]. When 
A(x) is evaluated using 


A (x) = dx' 


r j (xp -I 

Jx-x'L 


( 6 ) 


then selecting A x -0,A y = B ZQ x + A H corresponds physically to generating the uniform B zo by 
adding infinite current sheets located at* = ±°o, with current flowing in the ±y directions. In this 


case, only H and P y are independent constants of the motion. Use of this gauge is reasonable if 
B zo is produced primarily by distant cross-tail currents. Another simple gauge, in which A* - 
~ B z 0 y an ^ A y = A h , yields H and P x — ymv x + qA' x as constants of the motion. Physically, this 
suggests producing B zo by adding sheet currents at y = ±°o with current flowing in the Tx direc- 
tions. This gauge therefore is appropriate if B zo is produced primarily by currents flowing in the 
flanks. Since both forms of A yield identical magnetic fields, either can be used to follow particle 
motion. However, Chen and Palmadesso [1986] showed that even though P y and P' are both 
constants of the motion, they are not independent. Their Poisson bracket is nonzero \P' P 1 = 

y m ^zo’ ^ consequence of not having three independent constants of the motion is the possible 
appearance of stochastic particle orbits. 

Here we repeat the analysis of Chen and Palmadesso [1986] using B y (z) = B yo . The new 
vector potential can be written in simple gauges by adding either B yo z to A x or by adding -B yo x 
to A z . The former corresponds physically to the creation of B yo in the one-dimensional model by 
adding two current sheets located at z — ±°o with current flowing in the +x direction. This gauge 
suggests B yo is created by field-aligned currents in the mantle, or in the PSBL when only the 
inner CPS is modeled. The latter gauge corresponds to introducing sheets at x = ±°o with current 
flowing in the T z direction. Since there are two simple ways to introduce B zo by modifying the 
original Harris A and two equivalent ways to introduce B yo , A can be written in any of these 4 
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simple gauges. Two forms have A y = B zo x + A H , both of which yield P y = y mv y + qA y as a con- 
stant of the motion. The third form has A' x = B yo z - B zo y so that P\ = ymv x + qA\ is a constant 
of the motion. The fourth form of A yields no constant component of P. Since all 4 forms of A 
yield identical B fields, there are again a total of 3 constants: H, P y and P' Even though the 
above expressions for P y and P' x are different in the model with B yo = 0 than in the model with 
B yo * 0, the Poisson bracket [P' x , P y ] = ym& zo is the same. This suggests that stochastic particle 
motion is possible in both cases. We will see below that the most striking consequence of adding 
a modest B yo is a substantial change in the simple resonant behavior described in Section 1. 

E Fields and Acceleration 

For one-dimensional models with a uniform B z (z) - B zo , a uniform E y (z) = E yo in the 
Earth reference frame is transformed so that E' y = 0 as seen in a frame moving Earthward at u x = 
E yo /B zo , the deHoffman Teller frame. Our and many other analyses were done in this moving 
(primed) frame. When B yo = 0, all other components of E or B as measured in the Earth frame are 
equal (to order y- 1 ) to the components measured in the moving frame. 

When B yo * °> both E y an d E z differ significantly in the two reference frames. The moving 
frame was selected so that E' y = 0. Because guiding-center electrons and non-guiding-center ions 
are not distributed in the same manner along a field line, it was not possible to maintain charge 
neutrality throughout the current sheet unless a sheath-like electric field E\ was included [East- 
wood, 1974; paper 2]. With 3-keV protons and 3/8-keV electrons, the potential difference 
between z = 0 and z-L was found to be approximately 200 V in the moving frame for most cases 
studied. The resulting E\ ranges from approximately 0.03 mV/m for a quiet time current sheet 
with L = 1 R e to approximately 0.6 mV/m for an extremely thin current sheet with L = 0.05 R L 
In comparison, the steady E y0 in the Earth frame is approximately 0.15 - 0.3 mV/m if there is a 
uniform 50- 100 kV potential drop across a 50 R £ wide magnetotail. The Earth frame E z is given 
by Ez = y ( E 'z ~ u x B ' y ) ’ where y is nearly 1.0 for any reasonable magnetotail convection speed. 
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As a result, E z differs significantly from E z only when B y is nonzero. A relatively large u x is 
obtained by picking E yo = 0.3 mV/m and B zo = 3 nT, yielding u x = 100 km/s. If, for example, 
Byo ~ -1 nT, then the resulting difference between E z and E f z is = 0.1 mV/m. Even 

though the above figures show that E z can be substantially different in the two frames when B yo * 
0, the changes in E y and E z between the two frames leave E • B invariant. Both observers there- 
fore see the parallel acceleration of electrons that is required to maintain charge neutrality. 

Alekseyev and Kropotkin [1970] considered particle orbits near an infinitely thin current sheet, 
and Cowley [1978] used the constancy of P' x when B yo = 0, to show that the guiding center of a 
particle leaves the current sheet a distance 


vBJ 

A y = y 2 -y 1 = C| cos{X i| + | cosoc 2 R (7) 

zo 

in the cross-tail direction from the point the particle entered the current sheet. In (7), \B X \ and B 
are assumed to be the same at the entry and exit points, and a is the pitch angle. Repeating Cow- 
ley s derivation for the one-dimensional plasma sheet with B yo ^ 0 yields 


Ay = 


( V x2~ V x\) +Q yo( z 2~ z 0 


for the exact displacement of a single particle that enters the current sheet at z x with velocity v xl 
and leaves at z 2 , v x2 . Averaging over phase angle gives the mean (or guiding center) Ay 


Ay 


V r|^2 Cos a 2 | |^.rl COsa i ~ 

Kol + ~~ ~B[ 


+ ?r £ ( z 2“ z i) 


(9) 


Equation (9) reduces to (7) if the entry and exit points are at the same z. However, the last term in 
(9), which also can be written as B yo (z 2 - z { ) / B zo , is large if a particle enters in one hemisphere 
(e.g., zj > 0) and leaves in the other hemisphere (z 2 < 0). As noted above, field lines do not lie in 
the x-z plane when B yo ^ 0, so even a particle which exactly follows a field line enters and leaves 
the plasma sheet at substantially different y locations when it travels from one hemisphere to the 
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other. 

The addition of B yo and the resulting last term in (9) does not affect particle energization. 
Since E y = 0 in the moving frame, no energization is associated with motion in the y direction. In 
the Earth frame, E y * 0 so the extra Ay in (9) that results from the addition of B yo may appear to 
energize a particle that enters the current sheet from one hemisphere and leaves in the other. How- 
ever, the previously discussed extra E z term in the Earth reference frame (y 'ujt ), when multi- 
plied by (z 2 - Zj) , cancels the energization due to the extra Ay effect, at least to the order y- 1 
considered here. As a result, the addition of B yo does not affect particle energization in either the 
moving frame or the Earth frame. The above example points out the need to include both qAyE y 
and qAzE z when considering energization in the Earth frame when B yo * 0 rather than simply 
assuming that the acceleration of a plasma sheet particle is proportional to Ay. 

4. Orbital Effects of B y :N =1 Resonant Orbits 

The imposition of B yo modifies both the dynamics and the current-carrying properties of 

plasma sheet ions. Here these two effects are analyzed separately for ions at the k = 0.53 reso- 
nance. 

Average Pitch Angle Changes 

Figure 4 shows how the addition of B yo changes the dynamical nature of ion orbits near the N 
= 1 resonance. Between 100 and 300 groups of 100 protons each were followed to produce this 
and each of the other contour plots in this paper. Particles in any one group had the same energy W 
and moved in a fixed field model, so all particles in a single group had the same k. In Figure 4a 
(Aa), the rms pitch angle change, was evaluated using the untrapped ions in each group. All ions 
were injected from the northern hemisphere, so that the initial pitch angle oq > 90°. The quantity 
Aa was defined as |a f - a f | for ions that left the current sheet in the southern hemisphere so that 
the final pitch angle a f > 90°. The quantity Aa was defined as 1 180° - a. - a f | for ions that left in 
the northern hemisphere so that oy < 90°. With this definition Aa = 0° if an ion which began with 
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OC; 170 left with cx^ either 10 or 170°. This, and most of the following plots, are labeled by pro- 
ton energy in a given model field rather than by the K parameter so they could more easily be 
related to observations. Ions with W = 1, 2.8, and 5 keV have k = 0.69, 0.53, and 0.46 respec- 
tively. The other magnetic field parameters are fixed at B xo = 32 nT, B zo = 3.3 nT, and L = 1 .0 R E 
for all groups. Figure 4 a shows that (Aa) is less than 5° at the resonant energy of 2.8 keV when 
B yo = 0. One of the most striking observations of this work is that the region of small (Aa), 
which is characteristic of resonant orbits, completely vanishes by the time | £ | reaches 0.5 nT. 
The average pitch angle change is between 25° and 30° for all energies in this initially isotropic 
(T\\/T ± = 1) distribution when 0.5 nT < \B yo \ < 3 nT. In this range of |£ | , no evidence of com- 
pensated or resonant orbits appears at any of the energies shown in Figure 4. The second region of 
small (Aa), at large negative B yo , is discussed below. It should be noted that large (Aa) does not 
necessarily mean that there is a poor correlation between a^ and a,-. A large (Aa) often is pro- 
duced by a systematic difference between these two angles, as discussed in Section 6. 

Figure 4 b is a contour plot of R , the fraction of untrapped ions that both entered and left the 
plasma sheet in the northern hemisphere. Resonant 2.8 keV ions essentially all crossed the x = 0 
plane twice when B yo = 0. The ions appeared to simply bounce off or were all reflected by the 
current sheet, giving R = 1. When B yo was strong and negative the ions were guided directly 
through the current sheet with almost no net change in p or a (R = 0). A sample orbit illustrating 
this behavior is presented in the following section. 

^ ien Byo was strong and positive, ions injected at z > 0 usually mirrored several times before 
leaving the modeling region. Orbits were stochastic, with p and a changing erratically by a large 
amount each time the ion moved across the current sheet from north to south. Ions usually mir- 
rored in the southern hemisphere and then returned to cross the current sheet from south to north. 
For this south-to-north crossing in a field with large positive B yo the ions behaved the same as 
ions did when they crossed from north to south with a strong negative B yo , 


i.e. a and p were 
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almost unchanged. Ions therefore would mirror next in the northern hemisphere with almost the 
same trajectory they had just followed in the southern hemisphere. This sequence of mirroring at 
random locations in the southern hemisphere followed by an almost identical orbit in the northern 
hemisphere would continue until a particle happened to leave the southern side of the current 
sheet with a small enough pitch angle to escape. Figure 4b therefore shows that essentially all ions 
which escaped when B yo = +5 nT did so in the southern hemisphere (R = 0). 

Individual Orbits and {v y ) 

Rather than the many groups of several hundred protons each that were used to prepare Fig- 
ures 2 and 4, Figure 5 shows the orbit of a single 3-keV proton. One ion was injected at z = 1 R E 
with the same pitch and phase angles (a = 160°, \\f = 90°) into fields which differed only in the 
value of B yo . Other parameters were the same as in Figure 4. The radius of the spiral, and there- 
fore a and p, were almost the same when the ion entered and left the current sheet in the B =0 

yo 

panel. The final spiral radius, p, and a are substantially different in the B yo = -0.2 nT panel. The 
proton barely escaped after a single neutral sheet interaction, with a near 90°, when B yo - -0.4 
nT. The proton mirrored one or more times within the current sheet in each of the lower 3 panels. 
The resulting final cx is almost arbitrary in these latter cases, yielding a large Aoc. 

Even though p and oc changed substantially for \B yc ? | as small as several tenths of a nanotesla, 
the overall character of the orbit does not appear to have changed drastically in the top 3 panels. 
This suggests that the current-carrying properties, in which we are most interested, may not be 
quite so strongly dependent on B yo . 

Figure 6 shows the maximum (v y )/ v 0 , which usually is found near z = 0 (see Figure 2). This 
contour plot is based on the same orbits that were used to prepare Figure 4. The {v )/v 0 parame- 
ter still shows a weak increase at 2.8 keV when \B yo \ = 0.5 nT. The strong guiding action of a very 
large negative B yo again is evident, with (v y )/v 0 exceeding 0.8 at B yo = -5 nT. 

Figure 7 shows that the z dependence of (v y )/ v Q does not change drastically at the 2.8-keV 
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resonance energy when | B yo \ < 0.4 nT. This figure was prepared by tracing a group of 1000 pro- 
tons in each of the 6 field models with different B yo . Similar results were reported and discussed 
for the B yo = 0 case in paper 2, and by Ashour-Abdalla et al. [1992 b] using a two-dimensional 
magnetic field model. The magnitude of (v y )/ v 0 decreased, and the peak became slightly more 
confined near z = 0 as | B yo \ increased. The peak in j y (z) near but not exactly at z = 0, which is 
characteristic of resonant orbits, almost vanished when B y0 reached -1 nT. The j y (z) profile at 
Byo ~ ~ 1 in Figure 7 is more like a B yo — 0 profile dominated by stochastic orbits (k = 0.40, 

0.82) than by resonant orbits (Figure 2). 

5. Orbital Effects at Higher Resonances 

N = 2 Resonant Orbits 

Figure 8 shows that the situation is markedly different when W is increased to 23 keV, the N = 
2 or K — 0.32 resonance, in this same magnetic field. Even-A orbits involve an integer number of 
full sinusoidal oscillations during each current sheet interaction. An N — 2 ion that enters the cur- 
rent sheet from the northern outer CPS crosses the z = 0 plane 3 times, and leaves the current sheet 
to enter the southern outer CPS. The ions in Figure 8 came from an initial distribution with 
T \ / T ±. = 10 at the edge of the current sheet. Resonant effects in (Aa) were present, but much 
less dramatic, when a T { /T ± = 1 distribution was used. Figure 8 shows that there is a striking dif- 
ference between the effects of positive and negative B yo on even-A orbits [Dusenbery et al ., 
1993]. 

Figure 9 shows trajectories of two particles injected in the northern hemisphere. Ions gener- 
ally drift in the positive y direction in the current sheet. For large | B \ , it is possible to find the 
symmetric trajectories that are needed to yield small (Aa) only when B yo < 0. Ions injected from 
the southern hemisphere have small (Aa) when B yo is large and positive. Figure 9 also shows 
how large the cross-tail motion is in a small region near z = 0 when B yo is strong and negative. 
When B yo — —5 nT, even the N = 1 (k = 0.53) resonant particles, which were reflected from the 
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current sheet when B yo = 0, are guided directly through the current sheet (Figure 46). This is the 
region with veiy large (v y )/v Q in Figure 6. In the B yo = -5 nT field these k = 0.53 particles have 
trajectories that look very similar to the B yo = -2 nT orbit in Figure 9. A modest decrease in (Aa) 
also is evident near the point at which B yo = 3.7 nT in Figure 8. This structure is near the point at 
which B yo — B zo . However, an examination of individual orbits did not reveal any simple pattern 
such as those seen in the regions with very small (Aa). 

N = 3 Resonant Orbits 

Figure 1 0 shows the N = 3 or k = 0.23 resonance at W = 3.6 keV. This figure used B xo = 32 nT, 
B zo ~ 3-3 nT, L = 0.2 R E , and T^/Tj_ = 10. The region of very small (Aa) disappears at a few 
tenths of one nanotesla (Figure 10a). The N = 2 resonance is at W ~ 900 eV in this magnetic field. 
Figure 10a again illustrates that (Aa) is strongly dependent on the sign of B y at even-A reso- 
nances and nearly independent of the sign when N is odd. 

Figure 106 shows that the resonant structure is present, but relatively weak in its effect on 
cross-tail currents or (v y )/v 0 . This ratio exceeds 0.4 near B yo = 0 at both the N = 2 and N = 3 res- 
onances and drops below 0.2 half way between the resonances, near 2 keV. The guiding action of 
very strong B yo fields, illustrated in Figure 9, again produces much larger values of (v )/v 0 at 

Byo ~ than at B yo — 2 nT for these particles that were injected in the northern hemisphere. 

N = 6 and 7 Nonresonant Orbits 

Figure 1 1 shows that weak periodic structure can be detected near and even beyond k = K r . In 
Figure 11, which used B xo = 32 nT, B zo = 3.3 nT, and L = 0.04 R E , k decreases from 0.14 at W = 

1 keV to 0.088 at W = 6 keV, while K r = 0.10. Even though there are no regions of very small or 
very large (Aa), weak minima are evident when B yo = 0 at 1.8 keV (N = 6 or K = 0.12) and at 3.4 
keV (N - 7 or K = 0.10). Figure 11 was prepared using T^/T^ — 10 to emphasize the weak min- 
ima in (Aa), which are almost absent in a T^/T L = 1 distribution. It may be noted that the min- 
ima have shifted away from the locations predicted by (4) in this nonresonant region. The failure 
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of (4) to predict the exact locations of minimum (Aa) for nonresonant orbits is not surprising 
because these ions can bounce beyond the edge of the current sheet (Izl = L ) during the meander- 
ing portion of each orbit. The k parameter was defined [Buchner andZelenyi , 1989] using a linear 
approximation for bounce frequencies and radii of curvature near z = 0 rather than for particles 
that bounce out into a region of nearly uniform field. At veiy small k values, Buchner and Zelenyi 
[1989] showed that a normalized form of I z is more important than k in determining pitch angle 
diffusion rates. The previously discussed dependence of (Aa) on the sign of B yo for odd and 
even N also can be seen in Figure 11. A plot of (v y )/v Q showed essentially no dependence on W 

or k even when T { /T ± = 10. Cross-tail currents therefore show essentially no periodic indepen- 
dence when k < K r . 

We have used (3) to define k even when B yo * 0. If k 2 is defined as the ratio of the field line 
radius of curvature R c to the maximum particle gyroradius, both evaluated at z = 0, then a second 
parameter 



[Karimabadi et al ., 1990; Buchner and Zelenyi , 1991] also is important in determining orbital 
characteristics. The ratio Wk is 1.07, 1.26, 1.57, 1.97, and 2.45 when \B yo \ = 1, 2, 3, 4, and 5 
nT respectively, with B zo - 3.3 nT. Although phase space characteristics are dependent on or 
= K ( B yo /B zo ) [ Buchner and Zelenyi , 1991], Figures 4, 8, 10, and 1 1 show that the locations of 


resonances are nearly independent of B yo . When K > k,., the energy of an ion at the Mh resonance 
therefore is well predicted by (4), which involves k but not [Dusenbery et al ., 1993]. 


6. Low Altitude Signatures 


Figure 12 compares pitch angle changes produced by two processes. The basic parameter set 
used was B xo = 32 nT, B yo = 0, B zo = 3.3 nT, L = 1 .0 R E , and W = 2.8 keV. Either W or B yo was 
changed to prepare panels b, c, and d , as noted below. Each marker shows the initial and final 
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pitch angles of one ion that started and ended at Iz! = 4L/ 3 . The parameter ( 180° - a.) was plot- 
ted on the * axis because orbits started in the northern hemisphere, where incoming ions have oq > 
90 . The y axis shows the final pitch angle, o y. Figure 12a used the basic set of parameters, which 
placed ions near the N = 1 resonance (Figure 4). If compensation was exact for ions with all pitch 
and phase angles, then (180° — oc ; .) and oy would be equal. Figure 12a shows a narrow band of 
ions that very nearly satisfy this requirement for exact resonance. A second broader but still very 
sharply bounded distribution of final pitch angles is evident at large ( 180° - a.) . At these pitch 
angles, Aoc becomes phase angle dependent, with a well defined maximum Aoc . 

Figure 12 b shows the distribution of pitch angles for 3.5 keV ions, which are not near a reso- 
nance. Figure 4 shows that when B yo — 0, the mean pitch angle change for a 3.5 keV proton is 
approximately 15°. Figure 12 b presents more detailed pitch angle information for these ions. An 
initially field-aligned ion ( 180° - a. = 0°) leaves with Of = 15°. The final pitch angles are dis- 
tributed throughout a clearly bounded range for other oq. A simplistic picture that can explain the 
sharp cutoffs in oy for each oq (paper 1) is that the net effect of a current sheet interaction is to 
deflect all incident ion velocity vectors as a block through the same scatter angle experienced by a 
field-aligned ion (15° in this case). This model would confine points in Figure 12/? to a diagonal 
band with a height of 30° along the o tf axis. Figure 12 b shows that this simple model fails to 
explain two features of the calculations. First, for (180° — a.) > 30° the band of points exceeds 
the expected 30° range in o if, even though the band retains its sharp edges. This again suggests 
that Aoc has some phase angle dependence. In addition, the clustering of points in certain patterns 
within the allowed band is not explained. 

It is not surprising that Aa shows some dependence on (180° - a.) when this latter angle 
becomes large. Many of the analytic results have been derived for limiting cases, such as small I z 
and large or small k. One striking feature of orbit tracing analyses is that resonances such as that 
at k = 0.53 are so nearly independent of the injection pitch and phase angles [Buchner, 1991]. 
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Savenkov et al. [1991] and Brittnacher and Whipple [1991] studied the phase angle dependence of 
changes in I z . It was shown that I z of generally field-aligned ions jumps by a predictable amount 
when the ion crosses the first separatrix, from spiraling to meandering motion. The jump in pitch 
angle during the second separatrix crossing (from meandering back to spiraling motion) is 
strongly dependent of the phase angle of the incident ion. Changes in I z or pitch angle vary 
smoothly with phase angle, but at certain critical phase angles A I z becomes especially sensitive to 
the incident conditions. This behavior may explain the internal banding structure that is evident in 
Figures 12a and 12 b. 

The W = 5.5 keV protons in Figure 12c are farther from the resonance, so an initially field- 
aligned ion leaves the current sheet with c Xf = 52°. The pattern out to (180° -a.) = 10° is 
approximately described by a shift of all ion directions by 52°. In Figure 11c, some protons with 
(180° -a.) > 10° have Of very widely distributed. This behavior begins when the sharply 
bounded range of cy reaches 90°. At this point, ions begin to mirror within the current sheet rather 
than escaping after one current sheet encounter. These ions essentially are reinjected with much 
different pitch angles. Protons with oy> 90° in Figure 11c left the current sheet in the southern 
hemisphere. 

Figure 12 d was made by injecting the same 2.8-keV ions as in panel a into a field with B = 

L yo 

-0.2 nT. Figure 4a shows that both the (B yo =0,W= 3.5 keV) and the ( B yo = -0.2 nT, W =2.8 
keV) points lie on the (Aa) = 15° contour. A comparison of Figures 12 b and 12 d shows that mov- 
ing away from the ( B yo = 0, W = 2.8 keV) resonance point either by shifting W or by shifting B yo 
produces very similar pitch angle distributions near odd-A resonances. However, it should be 
recalled that odd-A resonances appear only when B yo = 0. It is not possible to return to an odd-A 
resonance by changing IT in a non-zero B yo . Buchner and Zelenyi [1991] previously compared 
changes in surface of section plots, and therefore in the type of orbit followed by an ion, that are 
produced by separately varying k and k s . This is equivalent to our changes of W and B . 
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It has been shown that low altitude satellites can detect ions in a broad energy range which 
precipitate from regions of the plasma sheet that are dominated by stochastic orbits [Ashour- 
Abdalla et al, 1992 a; Bosqued et al . , 1993]. A gap in which ion precipitation diminished also was 
seen at low altitude. The gap was attributed to a wall or Alfven-layer-like region at which k » 1 
that is associated with rapid drift across the tail (in a two-dimensional model) or around the Earth 
(in a three-dimensional model). Finally, it was suggested that groups of ions may be strongly 
accelerated near the neutral line to form beamlets, which also could be detected at low altitudes. 

We suggested in paper 1 that it may be possible to detect a different resonant orbit effect that 
would be observable only in very narrow energy bands. Some energy can be found at which parti- 
cles are resonant on almost any midtail field line. The requirement for the low altitude observation 
of a very narrow-banded feature is that field-aligned ions must suffer less than a 1 ° change in 
pitch angle as they pass through the current sheet. One consequence is that the loss cone will 
remain empty when equatorial pitch angles change by < 1°. Our finding that odd-A resonances 
completely disappear when a small B yo is added suggests that the proposed effect will not exist at 
odd-A resonances. However, very small pitch angle changes are possible at even-A resonances, so 
the proposed effect is investigated further here. 

Since we are interested only in ions that reach low altitude satellites, a highly field-aligned 
distribution (T l /T ± = 10 4 ) was injected at z = L to prepare Figure 13. The field parameters were 
^xo = 32 nT, B zo =3.3 nT, and L = 0.4 R £ . When B yo = 0, ions in an 0.1 keV wide energy band 
at the A = 2 resonance (3.66 keV) left the current sheet with less than 1° rms pitch angle changes. 
Figure 13 shows that a region with (Aa) < 1° remains when B yo becomes large. In fact, the 
energy band in which (Aa) < 1° becomes 2 or 3 times wider at the larger B yo values shown, sug- 
gesting that the effect might be easier to detect at low altitudes when B yo is large. Figure 13 again 
shows that the resonant energy depends only weakly upon B yo . 

If a resonant energy could be detected at low or middle altitudes, it would provide a way to 
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measure the ratio B 2 zo L/B xo appropriate to the Harris model at the equatorial end of the field line. 
Either of two specific effects may be observable. At any location, ions in the very narrow resonant 
energy bands would not be scattered enough to fill the usual upgoing loss cone. A downgoing loss 
cone therefore would exist, and could be seen at any middle or low altitude. Since downgoing ion 
fluxes generally are weak, it would be difficult to identify such a narrow band feature. However, 
intense upgoing ion beams are often observed. The lack of any pitch angle change in the current 
sheet would result in a very narrow energy downgoing ion beam. Ions at adjacent nonresonant 
energies would be scattered to form an isotropic rather than a beam-like downgoing distribution at 
low and middle altitudes. 


7. Birkeland Current Generation 

Overview 

This section uses the arguments proposed by Vasyliunas [1972] to connect cross-field drift and 
field-aligned current [Siscoe and Maynard, 1991; Heinemann and Pontius, 1991; Pontius , 1992; 
and Bit mingham, 1992]. The original analysis was based on the experimental observation that 
electrons carry nearly all the field-aligned current. This is true even within ion beams [e.g., Kauf- 
mann and Kintner, 1984]. As a consequence, the ion continuity equation 

V * (V B + \ di ) ] = -d (q t N.) Idt (11) 

requires that equal ion fluxes drift in and out the sides of a stationary flux tube during steady-state 
conditions. In (U),N i = N e = N is the number of ions and of electrons in a flux tube with unit cross 
section at the ionospheric end. The E x B drift velocity V# is the same for ions and electrons. 

Both V 5 and \ di , the sum of all other ion guiding center drift velocities, are referenced to the ion- 
osphere. 

A substantial net electron flux is permitted to drift through the sides of a stationary flux tube. 

Electron continuity is satisfied by electrons flowing into or out of the ionosphere through the end 
of the flux tube 



25 


V ’ [?A(V 5 + V^)] = -d(q e N e )/dt-j^ ( 12 ) 

where the electron q e is negative, and positive Birkeland current ;*n flows down to the ionosphere. 
We will consider only steady-state conditions in the remainder of this work, so d ( qN ) /dt = 0. 
Since the average ion is substantially more energetic than the average electron, \ de often is 
neglected relative to \ di and V# in the guiding center treatment. 

The above expressions show that the underlying cause of field-aligned current is that the 
derivatives of N e V de and A^V rf/ are different, but that N e = N t . In general, any difference between 
variations in fluxes resulting from ion and electron drift requires the presence of Birkeland current 
to maintain charge neutrality. We assume that any parallel electric fields that are needed to main- 
tain N e = Nj will form and either precipitate magnetospheric electrons or extract ionospheric elec- 
trons. 

In the present study, average drift velocities were calculated only inside the magnetotail cur- 
rent sheet. We therefore divide each flux tube into two segments. The current sheet segment is 
defined as extending between the z = 0 and z = L planes, where drift velocities have been calcu- 
lated. Beyond this point \ de - Y di = 0 in our one-dimensional model, though these drifts will be 
nonzero in a three-dimensional model. The flux tube is chosen to have a rectangular cross section 
at z — 0. The remainder of the flux tube extends from z = L to the ionosphere. This magnetospheric 
segment of the flux tube is selected to have a unit area at its ionospheric end. The guiding center 
approximations and the previous analyses are valid in the magnetospheric segment, except for the 
addition of current through the end at z = L. 

Current continuity and charge neutrality are maintained in the current sheet segment by con- 
sidering drift through the 4 sides of a flux tube and drift plus Birkeland current through the end at 
z = L. By symmetry, no current flows through the end at z = 0. Current flowing out the z = L end of 
the current sheet segment flows into the magnetospheric segment. Part of this current serves as an 
additional source to the usual guiding center analysis, and therefore adds to the current which 
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must reach the ionosphere. The analysis would be simpler if flux tube segments were defined with 
ends that are normal to B, but our calculations are all based on boxes aligned in the x, y, and z 
directions. It may be noted that E x B drift produces equal electron and ion fluxes, and no net cur- 
rent through the faces of the two segments at z = L. Neglecting the E x B drift contribution to cur- 
rents through these two z - L faces will not change the basic conclusions regarding Birkeland 
current flowing to the ionosphere. Since non-guiding-center aspects of particle motion are impor- 
tant only in the current sheet segment, this is the only segment that is considered in the remainder 
of this section. 

Elementaiy Models 

Figures 2 and 7, as well as several figures in paper 2, show that when E = 0 a resonant ion has 
a much larger (v^) than does a stochastic ion. This is true both for the peak (v ) and also for the 
average integrated through the current sheet. Electrons follow guiding center motion in this 
region. It is the difference between variations of the cross-tail drifts of guiding center electrons 
and of both resonant and stochastic ions that tends to create the charge imbalance that is studied 
here. The treatment in the present section differs from previous guiding center drift analyses 
through this emphasis on the abrupt changes in (v y ) that can be produced near resonances by 
modest changes in magnetic field parameters. 

Figure 14 illustrates the proposed mechanism. Both the original T89 model that was used in 
paper 1 and the version of T89 that was modified by Pulkkinen et al. [1992] in order to reproduce 
satellite observations during specific substorms showed that constant-K contours are curved in the 
sense shown by the solid curves near the Earth in Figure 14. The K labels used by Pulkkinen et al. 
are for electrons, but contours of ion k have the same shapes. The particular contours shown in 
Figure 14 a were selected because ions are stochastic at K = 0.82 and resonant at K = 0.53. 

The heavier and lighter arrows showing \ di at K = 0.53 and at k = 0.82 in Figure 14a are 
intended to indicate variations in cross-tail drift speed. Ions do not appear, disappear, or even 
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change their energy when viewed in the moving coordinate system of our one-dimensional model 
field, but simply change their average cross-tail drift velocities when they pass between resonant 
and stochastic regions. 

To provide the simplest model that can illustrate the effects of interest, we consider a region in 
which the magnetic field varies as one moves in the y and z directions, but is independent of x. We 
selected an electric field E y that is independent of y within the flux tube, but that may vary in x. 

The sketches all assume that V B is exactly in the x direction at z = 0 and that \ d is in the y direc- 
tion. 

Figure 1 5a shows particle fluxes flowing through the sides of a rectangular flux tube that 
extends from the ion k = 0.82 contour on the left to the K = 0.53 contour on the right. This orien- 
tation corresponds to the dusk (y > 0) side of Figure 14 a. Figure 15a is appropriate for a region 
containing a uniform E y and N that varies only in the y direction. Equal fluxes NV B of ions are 
shown drifting into the tailward side and out of the earthward side of the tube. Since V di is larger 
on the right side than on the left side of Figure 1 5a, the illustrated density gradient in the y direc- 
tion is required in order to satisfy ion continuity (11). 

Electrons drift earthward at the same speed and have the same density as ions, so the electron 
fluxes flowing into the tailward edge and out of the earthward edge also are equal in Figure 15a. 
Since the guiding center electron drift does not change rapidly in the y direction, more electrons 
will drift in from the side with higher density than drift out the side with lower density. An electric 
field must form so that the excess electrons will flow into the magnetospheric segment of the flux 
tube, contributing to an upward y,,. 

Current continuity also could be satisfied with the same B and perpendicular E fields used in 
Figure 15a by assuming that the density gradient is in the positive x direction. Ion and electron 
fluxes and drift velocities in this case are sketched in Figure 156. The uniform V B requires more 
electrons and ions to leave the earthward side of the flux tube where the density is high than enter 
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the tailward side. The ion continuity requirement (11) could sometimes be satisfied with no den- 
sity gradient in the y direction because ions drift more rapidly when they enter the dawn side of 
the flux tube than when they leave through the dusk side. Since V de is uniform, there is a net out- 
ward flux of electrons through the sides of the flux tube, which must be balanced by electrons that 
enter through the end at z — L. This contributes to a downward in the magnetospheric segment 
of the flux tube. 

As a final simple example, all the fluxes sketched in Figure 15 b can be obtained with N com- 
pletely uniform if E y increases as one moves in the earthward direction. An x-dependence of E y 
can produce the same change in foV B as does an x-dependence of N. There clearly are many com- 
binations of the elementary models sketched above that can satisfy (11) and (12) when V de and V di 
have different y-dependence. Almost all such models require Birkeland currents. 

Figures 15c and 1 5d illustrate the same simple models applied to the dawn side of Figure 14a. 
Figures 15a and 15c both have the cross-tail and earthward ion drift currents separately satisfying 
the steady-state continuity equation. All Birkeland current is generated by the change in cross-tail 
electron flux that is associated with the density and accompanying pressure gradient. Since V de is 
relatively small, this model results in a modest field-aligned current density. It may be noted that 
Figure 5 of Kaufmann et al. [1993/?], which is based on the T89 model, shows little y dependence 
of the cross-tail current after integrating through the thickness of the current sheet. 

In contrast, cross-tail and earthward ion currents each vary in Figures 15 b and \5d. Equation 
(1 1) is satisfied because the two components of the divergence cancel. In this case Birkeland cur- 
rent arises because of the change in NV B , the earthward electron drift flux. This change in NV B can 
be produced either by an x-dependence in E y or a gradient of N in the x direction. Since V di is 
large, the change in NV di and therefore of NV B and /|, in Figures 1 5b and 15 d may be large. Note 
that Figures 15 b and 1 5d require different gradients in N or E y on the dawn and dusk sides. The 
product NV b must increase earthward on the dusk side and must simultaneously increase tailward 
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on the dawn side of the magnetotail. As one example, this situation is possible even if N is uni- 
form and the total cross-tail potential drop is independent of x. In this case, the region of strongest 
E y would have to shift from the dusk side to the dawn side as one moved out the tail to be consis- 
tent with Figures 156 and 15 d. 

Current Magnitudes 

The magnitudes of the fluxes described above can be estimated to see if it is reasonable to 
drive a significant night side steady-state Birkeland current solely from the magnetotail current 
sheet. Iijima and Potemra [1978] found that the average ionospheric yj| is approximately 1 
^ A /m . With B = 50,000 nT in the ionosphere and B zo = 5 nT in the equatorial plane, a (not nec- 
essarily rectangular) 1 m 2 area in the ionosphere projects to a square box, similar to those 
sketched in Figure 15, with a = 100 m sides. Assuming an average number density of n = 0.3 

_3 . 

cm in the central plasma sheet [Baumjohann et al., 1989] and a current sheet with half-thick- 
ness of h = 1 R e yields a flux of nha{v di ) = 2 x 10 14 (v di ) ions/s passing through the left and 
right sides of the flux tube. The quantity (v di ) is v y (z) averaged through the current sheet for the 
geometry in Figure 15. Using the model in Figure 156 as an example, the rate at which ions accu- 
mulate in the flux tube owing to cross-tail drift is 2 x 10 14 a (d(v di )/dy) = 2 x 10 16 ( d(v di )/dy ) 
ions/s • As discussed above, Figure 156 shows a net loss of both ions and electrons through the 
other two sides at this same rate owing to E x B drift. The net effect is for the electron content of 
the flux tube to decrease at a rate of 2 x 10 16 <7 (v di )f dy electrons/s . The corresponding expression 
for the model in Figure 15a can be written as 2 x 10 16 ((v de )/(v di )) d(v di )/dy electrons/s. An 
ionospheric y N of 1 jiA /m 2 requires 6 x 10 12 electrons/s flowing up from the ionosphere through 
the end of the flux tube. The value of d(v di )/ dy required in Figure 156 for steady state therefore 
is 3 x 10" 4 (m/s)/m or 2 (km/s )/R £ . In Figure 15a d(v di )/dy must be ( {v di )/{v de )) times larger. 

To see if the above d{v di )/dy estimates are reasonable, note that Ampere’s law requires a 
total cross-tail current, integrated through the current sheet thickness, of K y - 32 mA/m if B x is 
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20 nT on one side of the current sheet and —20 nT on the other side. Assuming that most cross-tail 
current is carried by ions, this K y * 2 enh(v di ) requires (v di ) « 50 km/s using the same density 
and half-thickness as before. The above estimates therefore show that (v di ) in Figure 156 must 
change by several percent per Earth radius of displacement in the y direction to supply all the 
electrons needed for a 1 pA /m 2 ionospheric Birkeland current. 

Drift Paths 

We traced several ion groups in the full three-dimensional T89 Kp = 4 magnetosphere model. 
The purpose was to see if (v di ) • V k is large enough so that most ions would drift from regions 
with k = 0.40 to 0.53 or from 0.53 to 0.82 before reaching the flanks. We chose E = 0 to isolate 
cross-tail drift from ExB drift effects, and from the changes in K due to particle energization. 
The use of E = 0 also simplifies the analysis by making it unnecessary to evaluate ExB drift 
through the boundary between the current sheet and magnetospheric segments of a flux tube. Note 
that the model in Figures 15a and 15c can be drawn with E = 0 and therefore NV B = 0. 

A group of 5 keV protons was injected at x = -20 R E , y = z = 0 in the T89 model to produce 
Figure 16. The k parameter was calculated each time an ion crossed the equator. The k parameter 
was defined as (R c / p z ) 1/2 where R c , the field line radius of curvature, was calculated using 
= - (B • V) B / B 2 at each equatorial plane crossing and p z = mv/eB z was evaluated at z 
~ 0. Since E = 0, W and v are constant, and changes in k are produced by drift to regions with dif- 
ferent equatorial B z and different R c . 

Figure \6a is a plot of equatorial crossing points, and shows that (v di ) was primarily in the 
cross-tail direction for lyl < 10 R E . Drift paths became more nearly circular, centered at the Earth, 
beyond lyl = 10 R E , so a comparison with one-dimensional results becomes less meaningful here. 
Figure 166 shows how k changed as the ions moved across the tail. This group of ions moved 
approximately 3.5 R E in the y direction to go from k = 0.53 to K = 0.82. These results are similar 
to those shown in Pulkkinen et al. [1992] and paper 1. The d{v di )/dy of 2 (km/s )/R E needed to 
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generate j n = 1 pA /m 2 therefore requires a 15% change in (v di ) between the k = 0.53 and k = 
0.82 contours. A change of this magnitude appears reasonable from our one-dimensional analysis 
(e.g. Figure 2), but a full three-dimensional calculation will be required before definitive conclu- 
sions can be drawn. At present, the suggestion that significant field-aligned currents are associated 
with cross-tail drift between regions dominated by stochastic orbits and regions dominated by res- 
onant or guiding center orbits must be considered speculative. 

8. Conclusions 


Odd-N Resonances 


It was shown that the addition of a uniform cross-tail magnetic field, which is comparable to 
that typically measured in the magnetotail, completely eliminated the resonances that exist when 
B yo = 0 for particles that left the current sheet on the same side from which they entered. Pitch 
angle changes were found to depend especially sensitively on B yo near these odd-A resonances. 
As a result, all ions that flow up from the ionosphere with energies that place them either near an 
odd-TV resonance or in a stochastic region between resonances will suffer large pitch angle 
changes in the typical current sheet. Upgoing beams of all such ions will be so strongly deflected 
at the equator that they will not appear beam-like after a single current sheet interaction. Simi- 
larly, an empty loss cone in an upgoing distribution of such ions will be filled during each interac- 
tion with the current sheet. The electric current carried by odd-A resonant ions was found to 
depend a little less sensitively upon B yo than did the average pitch angle change. 

Even-N Resonances 


Ions were found to behave differently near even-A resonances. Orbits can remain more sym- 
metric at these energies (Figure 9 b\ so that resonant behavior is not so strongly affected by the 
addition or removal of a typical B yo (Figures 8 and 10). Narrow energy bands were found in 
which only very small pitch angle changes took place during an even-A resonant current sheet 
interaction (Figure 13). Ions near these resonances suffer little scattering when they enter the cur- 
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rent sheet from the northern hemisphere and leave in the southern hemisphere if B < 0. Particles 
going the other way through the current sheet suffer little scattering when B yo > 0. Unscattered 
ions can be recognized at low or middle altitudes by the presence of downgoing loss cones (loss 
cones that have not been filled in during traversal of the current sheet) or by the presence of down- 
going ion beams (upgoing beams which have passed through the current sheet). These features 
will be restricted to very narrow energy ranges. Detection of either signature at low or middle alti- 
tudes should provide information about the equatorial magnetic field structure. Since cross-tail 
current depends less strongly upon B yo than does pitch angle scattering, the structure of the cur- 
rent sheet appears to be changed only modestly by the addition of a typical B yo near an even-iV 
resonance. 

Cross-Tail and Field-Aligned Currents 

The study of currents carried by individual charged particles provided new insight into possi- 
ble causes of some steady-state Birkeland currents. The feature that is required for operation of 
the proposed mechanism to drive field-aligned currents is for cross-tail drift to carry ions across 
contours of constant k. When (\ di ) • V k is nonzero in the near-Earth plasma sheet, i.e. when con- 
tours of constant k are curved as in Figure 14, then ions drift between regions characterized by 
different dominant orbit types. The lowest altitude resonant and stochastic regions are most 
important because there is a large energy separation between ions at adjacent resonances (Figures 
2, 4, and 8). The large energy separation permits most ions in a Maxwellian distribution to follow 
resonant orbits at one location in the tail, and to follow stochastic orbits at a different y. 

One-dimensional orbit calculations showed substantial changes in the average drift speed dur- 
ing the transition between points characterized by k = 0.53 and by k = 0.82 when B yo = 0. When 
B yo is substantial, the first resonant transition is that between k = 0.32 and K = 0.40. Similar 
changes also exist at the k > 1.5 guiding center to k = 0.82 stochastic transition. These are the 
points at which integral ion densities and pressures may change, and at which Birkeland currents 
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may be needed to maintain charge neutrality. We have not carried out the three-dimensional calcu- 
lations that would be needed to compare the proposed mechanism to results from the usual guid- 
ing center analysis. 

The calculations assumed that the same flux of ions drifts into and out of the sides of a station- 
ary magnetic flux tube. This assumption is commonly used because field-aligned currents are 
observed to be carried predominantly by electrons. Ions do not enter or leave the plasma sheet, but 
merely drift faster or slower as K changes. The other standard assumption is that the plasma is 
neutral. Since electrons follow guiding center motion in this magnetotail region, their drift veloc- 
ity does not change abruptly near an ion resonance. As a result, the electron flux entering the sides 
of the magnetic flux tube does not equal the flux leaving the tube. The difference must be supplied 
by electrons entering from or precipitating into the ionosphere through the tube ends. The possi- 
bilities that electrons accumulate primarily through changes in the cross-tail fluxes (Figures 15a 
and 15c) and through changes in the flux associated with earthward E x B drift (Figures 1 5b and 
\5d) both were considered. In either case, electrons carry the current that is diverted to the iono- 
sphere, even though ions carry most of the cross-tail current. Since a specific flux of electrons is 
required to maintain neutrality, the proposed mechanism operates as a current generator rather 
than a voltage generator. Parallel voltage drops must form, but their magnitudes will be whatever 
is needed to drive the required y ( ( . 

The present analysis does not address the question of why constant k surfaces are curved as 
they are in the standard magnetic field models, and only steady-state conditions were investigated 
in this work. The ionosphere and its response to imposed currents are critical, especially when 
temporal effects are considered. These essential but very complex aspects of the Birkeland current 
system have not been addressed in the present work. 

It should be noted that the abrupt changes in (v di ) predicted near ion resonances by orbit trac- 
ing studies will not appear in simulations that use MHD or the guiding center equations. The pro- 
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posed Birkeland currents should appear in any kinetic simulation that bases electric currents on 
the actual motion of individual particles, even if orbit tracing is relatively crude. In this regard, the 
present analysis may be compared with the interesting conclusion of Wolf and Pontius [1993] that 
in certain circumstances the drift of non-guiding-center particles, when averaged throughout a 
flux tube, is the same as the drift that would be calculated using the guiding center equations. 
Their analysis is based on the assumption that the distribution function can be written as a func- 
tion of only W and P y , which are constants of the motion. The total energy W is constant in the 
steady state analyses that we consider, and it was noted in Section 3 that P y is constant for one- 
and two-dimensional analyses. However, a requirement for variations of B or E in the y direction 
is inherent in any three-dimensional model that includes Birkeland currents associated with cross- 
tail drift resonances. As a result, P y will not be a constant of the motion. Our one-dimensional 
analysis (paper 2) showed that the required self-consistent fly) is strongly dependent on variables 
other than W and P y or v y . More realistic two-dimensional analyses show a substantial x, z off- 
diagonal component of the pressure tensor [Burkhart et al , 1992]. Further work clearly is needed 
before definitive conclusions can be drawn regarding the possible association between non-guid- 
ing-center motion and Birkeland current. 
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Figure Captions 

Fig. 1 . Particle orbits that follow invariant tori (a) ring-like torus, (b) half of a split torus, (c) full 
split torus. 

Fig. 2. Dependence upon k of the normalized cross-tail drift velocity distribution (v^ (z) )/v 0 
when B yo - 0. 

Fig. 3. Sketches of B y produced by different sources. 

Fig. 4. {a) Root mean square pitch angle changes for protons resulting from current sheet interac- 
tions near the iV = 1 resonance (b) Fraction of the protons that left the current sheet on the same 
side from which they entered. 

Fig. 5. Sample proton orbits in fields with different B yo . 

Fig. 6. Peak value of the normalized cross-tail drift velocity near the TV = 1 resonance. 

Fig. 7. Normalized distribution of the cross-tail drift velocity for several B yo . The ion energy cor- 
responds to the TV = 1 resonance. 

Fig 8. Root mean square pitch angle changes for protons resulting from current sheet interactions 
near the TV = 2 resonance. 

Fig. 9. Sample proton orbits in model fields with B yo = ±2 nT. 

Fig. 10. {a) Root mean square pitch angle changes for protons resulting from current sheet inter- 
actions near the TV = 2 (W = 900 eV) and TV = 3 (W = 3.6 keV) resonances, (b) peak value of the 
normalized cross-tail drift velocity. 

Fig. 1 1 . Root mean square pitch angle changes for protons resulting from current sheet interac- 
tions near the TV = 6 (W = 1 .8 keV) and TV = 7 (W = 3.4 keV) resonances. 

Fig. 12. Pitch angles for individual protons before (oq) and after (c y) a current sheet interaction. 
Fig. 13. Root mean square pitch angle changes for protons resulting from current sheet interac- 
tions near the TV = 2 resonance. 
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Fig. 14. Sketches of contours of constant k, the ion and electron cross-tail drift velocities V di and 
V de , and the number of ions and electrons in a unit flux tube N. 

Fig. 15. Average cross-tail drift velocities V di and V de , flux tube content N, and fluxes of ions and 

electrons for simple models of the transition region between stochastic and resonant orbits, as 
shown in Figure 14 a. 

Fig. 1 6. (a) x, y locations of z = 0 plane crossings for a group of ions in the three-dimensional T89 
K P = 4 model (b) The k parameter evaluated at each z = 0 plane crossing. 
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